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ABSTRACT 

We study a quantum particle propagating through a “quantum me¬ 
chanically chaotic” background, described by parametric random matri¬ 
ces with only short range spatial correlations. The particle is found to 
exhibit turbulent-like diffusion under very general situations, without 
the apriori introduction of power law noise or scaling in the background 
properties. 
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The diffusion of particles in a complex background, such as disordered media, has often an 
anomalous or Levy type character . In contrast to classical diffusion, where the variance 
of the displacement of a particle or tracer grows linearly in time, anomalous transport can 
have ( R 2 ) oc t a , with q^I, although it is not limited to power law behavior. The diffusive 
properties of such systems can hence range from enhanced to dispersive dynamics. The 
microscopic understanding of such stochastic processes in terms of the underlying (chaotic) 
dynamics remains an active open area of study |2|. While the conventional limits of Levy 
processes are 0 < a < 2, there exist superdiffusive dynamics beyond this range. One of the 
extreme examples is turbulent diffusion, which occurs when the background is turbulent, 
originally put in evidence in atmospheric measurements [|3]. In such situations, the average 
separation R between two tracers can increase as fast as (R 2 ) ~ t 3 . Richardson [[| postulated 
a Fokker-Planck equation of the form: 


dV(R, t) 
dt 


_d_ 

dR 


V(R) 


&P(R,t) 

dR 


V(R) oc R 4/3 , 


( 1 ) 


which by design reproduced his measurements. Kolmogorov, in studies of turbulence, pro¬ 
posed an energy-wavenumber scaling law (known as the |-law) and further suggested that 
dissipative behavior is spatially dependent [|. Refinements of these scaling arguments in¬ 
corporated intermittent corrections into the power law behavior [|(]. More recently, a new 
class of random walks, termed Levy walks, have incorporated Kolmogorov’s scaling to derive 
Richardson’s t 3 law [0 (including intermittency corrections), showing further that the scal¬ 
ing does not necessarily imply the latter. One of the common assumptions in the description 
of anomalous transport is a power law behavior of some input distribution function. For 
instance in Langevin or Fokker-Planck approaches, a power law behavior is generally chosen 
for the distribution of thermal noise [7j or in the spatial correlations [|lj; in deterministic 
chaotic models, power-law amplification is used in the random walks, long algebraic 
tails in step distributions or sticking times are used in fractional diffusion equations, 
the power of the fractional derivatives are chosen to describe the anticipated behavior |5]. 
This is not to say that such power-law behavior is not seen. Indeed experimentally one can 
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justify some of these assumptions ['ldj . 

It is instructive to preface our analysis with the results from the classical diffusion prob¬ 
lem in d— dimensions in the presence of a quenched random force F fl {x) and thermal activa¬ 
tion rj fli (t). Consider the Langevin equation [[TJ 

dx 


dt 


= F^x) +77 /1 (t), 


( 2 ) 


where spatial averages are F^ix) = F 0i/i , [F^x) — Fq^[F v {x') — F 0>v ,\ = G^x — x') , and 


the time average of the noise term is 7^ (£)?"/„ (t') = 2 D5(t — t')8^ v . One of the crucial issues 
is the construction of the statistical correlation function G^(x). By using a power law 
dependence for this function, it can be shown that the dynamics displays diffusion in three 
dimensions. Thus, the long distance correlations, which are characteristic to superdiffusive 
dynamics, have to be incorporated into the problem form the very beginning [|l]. 

In this letter we would like to approach the problem from a slightly different perspective 
and extend this type of study to the quantum regime, using the Schrodinger equation to 
describe the dynamics, instead of the Langevin formulation with the thermal activation. The 
fluctuations will emerge from the dynamics of a test particle in the presence of a correlated 
chaotic quantum background. We will see that turbulent-like behavior can be manifest on 
certain time scales under fairly general conditions and that power-law assumptions are not 
necessary for this. 

Since it is known that chaotic dynamics can induce superdiffusive behavior in classical 
Id-maps ||, as well as in classical motion coupled to quantum backgrounds JTl|, a natural 
starting point is to utilize random matrix Hamiltonians, which are essentially the quantum 
counterparts of classical chaotic systems. This will provide both a reasonable physical picture 
as well as a tractable framework for the analysis of the diffusion process. We take a model 
space which is the direct product of the Hilbert space of the test particle with position R, 
and the finite dimensional (albeit large) background space defined by a complete basis of 
states |i), with i = 1 1. The background, denoted Vij(R) , will be taken to be 

quenched (time-independent), and chaotic in the sense that the spatial inhomogeneity is 
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described by a deformed, parametric, banded random matrix: 


Vij(R) — Uo t ij + Ui t ij(R) 


( 3 ) 


where U\^j{R) is a real symmetric matrix, and an element of the Gaussian orthogonal 
ensemble (GOE) |12j]. As such, it is characterized by the first two cumulants: 


Ui,ij(R) — 0 


( 4 ) 


Ui,ij(R)Ui,ki(R') — [SikSji + SuSjk]Gij(R, R'). 


Here the overline indicates the average over the GOE |L3|. The density of states of the 
background is defined through the diagonal matrix Uo : ij = with the constant average 

level density po (po(^i+i — G*) = 1). The correlation function Gij(R,R’) is parameterized as 


pi 

Gij (-R, R') = -—exp 

lirpo 


(fy ~ GLjf 
2 «§ 


G 


'R-R ,s 

Xn 


( 5 ) 


This incorporates a bandedness for the random matrix, with an effective width N 0 = 
serving to limit the interaction range to the nearby states, an overall strength denoted T^ 
(also known as the spreading width [[[I|), and a spatial correlation function G(R/X 0 ) with 
length scale A"o, normalized such that G( 0) = 1. We will assume that the statistics of the 
background are translationally invariant, Gij(R,R') = Gij(R — R'), although this is not 
crucial. This function is the matrix analog of the correlated noise used in Eq. (2). But 
instead of building in long power law tails, we will use G(x) = exp[—x 2 /2], which provides 
a rapid spatial decorrelation. Next we couple a test particle of mass M to this “chaotic” 
background through the Hamiltonian: 

n 2 




( 6 ) 


Here, as in all formulas we present, R can be interpreted as a variable of an arbitrary 
dimensionality, even though some of the formulas we write explicitly for the 1-d case. Using 
the Feynman and Vernon |TB| formalism one can represent the density matrix for our test 
particle through the following path integral formula 


4 














rR'(t)=R' 


x exp ^ — [S 0 (R(t)) - S 0 (R’(t ))] j. F(R(t), 


p(R,R',t ) = / dX 0 dY 0 p 0 (R 0 ,R' 0 ) / / 2?i2'(t) 

J Jr(o)=Rq Jr'(o)=r ' 0 


(7) 


where So(-R(f)) is just the classical action for a free particle and po{Ro, R' () ) is the initial 
density matrix of the test particle. Taking advantage of the large Nq limit, it is possible to 
explicitly compute the influence functional for our Hamiltonian in the adiabatic limit |T7|,T5| 
(which we discuss below), where we find 


F(R,R',t ) = exp j\lt' [G([R(t') - R'(t')\/X 0 ) - 1]|. 


( 8 ) 


Thus the density matrix for the test particle satisfies the following equation |18 

n 2 


ih^p(R,R’,t) = 


2 M 


(d 2 R -d l2 ) + iVt(G(R,R')-l)) 


p{R, 


(9) 


In terms of the new variables r = (R + R')/ 2 and s = R — R' the solution of this equation 
can be found through quadratures (easily verified by direct substitution): 


p(r, s,t) — / dr' 


dk 


kt 


M P ° r ’ S “M exp 


ik{r — r') f s 


h 


hk Js-kt/M 


ds'[G(s'/X o) - 1] 


( 10 ) 


where P o(r,s) = p(r,s,t = 0) is the initial density matrix at t — 0. If we take the initial 
state to be a Gaussian 'ipo(R) = exp[—i? 2 /4a 2 ]/[27rcr 2 ] 1,/4 , the initial density matrix is 


p a (R, R’) = -X=e- = _J =e -(4r> + ,>)/8^ 


V2 


7TCT 


V 27TCT 2 


( 11 ) 


The adiabatic condition in which the influence functional (8) is valid, restricts the velocity 
V of the test particle such that the time scale (. Xq/V) is no greater than that of the back¬ 
ground characterized by h/n 0 : V max ~ n 0 X 0 /h. This can be used to constrain the average 
momentum of the initial wavepacket, through the width a. For our initial gaussian, these 
are related by ( P 2 ) = h 2 /(4cr 2 ) ~ (. MV ) 2 , so we require : 

h 2 


G > G m m ~ 


2,Ad XqRq 


( 12 ) 
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To extract the diffusive properties of the wavepacket, we compute the cumulants of the 
coordinate R directly from the coordinate distribution V(R,t ) = p(r, s = 0 ,t). This is done 
by constructing the characteristic function for the coordinate distribution, d(k,t), defined 
by taking the Fourier transform of V(R,t): 


d(k,t) = I drpo(r,—kt/M) exp 


ikr Y^M 


H + hk 


kt r° 

M J-kt/M 


dsG(s) 


(13) 


= exp <-- 


1 ( a k' 

n 


k 


2 V 2 Mo 


' t2| r+MJW2 ” (-1)" f 

kh ^ 


kt 


^ n!(2n + 1) \^2MX 0 ’) 


2n+l' 


= exp 


£ < ( Rm )> 


. 771—1 


ml 


All the cumulants are easily identified. The second cumulant, which measures the spreading 
of the wavepacket, is given by 


R" 


r ol 2 

= J dRR 2 V(R,t) = -h 2 —d{k,t ) 


k =0 


= R + 


n 1 


7 1 2 + 


r l h 




(14) 

(15) 


4M 2 ff r ' 3 M 2 Xl 

The terms are readily identified. The first is the initial width at t — 0, and the second is the 
natural spreading of Eq. (pT|) due to free expansion, which is the only dynamical contribution 
when the background is removed (E*- = 0). The dissipative contribution which arises from 
the background displays the diffusion associated with turbulent backgrounds, namely the t 3 
character. One can see that the turbulent-like contribution becomes dominant on the time 
scale 


t T « 


3 x 2 n 

4a 2 T! ' 


(16) 


The momentum distribution V(P,t) and its characteristic function D(s,t), are given by 


V{P,t) = 

JdRdR'ex p( lP ^ R f R 

-'J p(R, R',t) 

(17) 

= 

r , {IPs' 

J ds exp (a , 

) D(s,t), 


(18) 

<E > - 


O 

II 

CO 

nr 1 

“1“ 0 t. 

x 2 0 

(19) 
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One can see from Eq. flT9|) that in the absence of coupling to the background (F* = 0), the 
momentum cumulant is constant and given by the usual value for a wavepacket. The coupling 
to the background makes the momentum variance increase linearly with time. Because this 
turbulent-like behavior is limited to the adiabatic regime, the maximum time scale for 
turbulent-like diffusion to be present is given by the condition ( ( P 2 ) ) 1 ^ 2 
MX 0 K 0 /h, or: 


MV max ~ 


t 


max ~ g f T 


cr 


cr n 


( 20 ) 


which depends only upon the initial width of the wavepacket a. Hence for times on the 
scale t? < t < t max , the diffusion of the wavepacket will have a turbulent-like character. 
(i m ax > £t requires only that a > 2 cr min ) . For t > f max , the character of the interaction with 
the background changes over to a diabatic behavior, where the above form for the influence 
functional is no longer valid. This is not to say however that the dynamics ceases to have a 
turbulent-like behavior, only that our adiabatic expression for the influence functional (8) 
has a limited range of applicability. 

Similar to the Langevin approach in (2), the anomalous diffusion arises from the proper¬ 
ties of the spatial correlations, but for quite different reasons. Because we are using a random 
matrix ensemble to model the background properties, the correlation function G{x) must be 
positive definite, and can be classified by its short distance behavior: G(x) Ril* c\x\ a + • • • 
where the range is restricted to 0 < a < 2 and x = R/X o |I3|. If we want a smoothly 
correlated background then a = 2, while for a Brownian motion type spatial fluctuations 
a = 1. We will only consider here the case of smooth spatial correlations, a = 2. (The case 
a < 2 would be interesting to consider further, as the diffusion would be characterized by 
very long spatial tails and infinite moments, analogous to the Levy stable laws |lj.) The 
Gaussian correlation function provides generic results for any a = 2 correlator, which can be 
seen from the definition of the second cumulant. As the inverse is also a Gaussian, neither 
the spatial correlation nor its Fourier transform exhibit any long range correlations. All 
long range diffusive behavior emerges from spatial inhomogeneities on the scale Xq. If we 
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vary this scale, taking the limit X 0 —> 0, the anomalous diffusion is enhanced, as can be 
seen directly from ( ( R 1 2 ) ). The opposite limit, A" 0 —*► oo, corresponds to a constant random 
background (£/i,p(-R) is replaced with a fixed random matrix), and the turbulent diffusion 
vanishes, since the spatial domain on which it is active is never reached. 

By examining the dynamics of a wave packet in a chaotic background, we have found that 
turbulent-like diffusion can emerge under very general circumstances, with only the input 
of the short distance spatial correlations in the background on a finite scale Xy. no power 
law distributions are assumed. Further, our results for the full coordinate and momentum 
distributions, V(R,t) and V(P,t), do not exhibit the usual scaling behaviors or power law 
properties used in previous studies. Our dynamics is a statistical limit which emerges from 
the random matrix solution to the influence functional, and as the fluctuations are gaussian, 
cumulants higher than second order are not invoked. As with classical turbulent diffusion, 
which can be generalized to included intermittency corrections and so forth, a more general 
class of this turbulent-like quantum diffusion can be explored by considering various types of 
backgrounds. This would include, for example, stochastic rather smooth spatial correlation 
functions G(x) characterized by a < 2, corrections to the density of states for non-constant 
behaviors, or inclusion of higher cumulants in the background. In addition, the role of Ti, 
in particular, how the turbulent diffusion survives the limit h —> 0. These might provide 
a more general formulation of the quantum analog of diffusion in chaotic backgrounds, in 
which a classical limit might eventually recover intermittency corrections. 
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